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Abstract We study the interplay between disorder, interactions and 
decoherence induced by spontaneous emission process. Interactions are 
included in the Anderson model via a mean-field approximation, and a 
simple model for spontaneous emission is introduced. Numerical simu- 
lations allow us to study the effects of decoherence on different dynam- 
ical regimes. Physical pictures for the mechanisms at play are discussed 
' and provide simple interpretations. Finally, we discuss the validity of 

scaling laws on the initial state width. 



1 Introduction 

The Anderson model [T], is a simple, tractable, model describing disorder effects on 
electrons in a crystal. The model predicts that in disordered crystals the eigenfunc- 
tions may become spatially localised Q, in sharp contrast with the Bloch functions of a 
perfect crystal, a phenomenon called Anderson localisation. However, due to its very 
simplicity, the model (in its original form) overlooks potentially important effects: It 
is a single particle, zero temperature model. 

Recently, analogs of the Anderson model have been realized with cold or ultracold 
atoms, which can be made much closer to the Anderson model than a crystal. A direct 
translation of the Anderson model can be obtained by placing ultracold atoms in a 
far-detuned laser speckle, which realizes a random mechanical potential affecting the 
center-of-mass of the atom, and allowed observation of the localisation in one dimen- 
sion (ID) with bosons |2I3| . The 3D localisation was also studied with both fermions 
[1] and bosons [5] . The quasiperiodic kicked rotor can also realize an equivalent of both 
ID |6I7I8) and 3D systems |9I10| : the latter allowing a detailed experimental study 
of the Anderson transition, including measurements of the critical wavefunction 
and of its critical exponent [10112113] . Ultracold atom systems allow a large control 
of decoherence, which is essentially due to spontaneous emission, and of atom-atom 
interactions, for which an experimental "knob" is provided by the Feshbach reson- 
ances [14], making them an ideal tool for studies of the interplay of quantum effects, 
disorder, interactions, and decoherence. 



a http://www.phlam.univ-lillel.fr/atfr/cq 

1 In ID all eigenfunctions are localised, in 3D an energy "mobility edge" separates localised 
from delocalised ones. 
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We have recently studied numerically and theoretically the influence of interac- 
tions in a bosonic ID Anderson system |15j . which is a very active research domain 
|16I17) . The sensitivity of the system to the initial state was shown to be character- 
ised by a scaling law depending on the initial state. In the present paper we introduce 
the effect of decoherence due to spontaneous emission in the problem, and show that 
Anderson localisation (AL) is promptly destroyed by it, but the resulting dynamics 
still respects the scaling law mentioned above. This confirms the pertinence of the 
scaling strategy, and indicates that it is a relevant approach for studying disordered 
interacting systems. 



2 The model 

We essentially use in the present paper the same model as in |15j . except for the 
introduction of decoherence effects as discussed below. Briefly, we start with a gas 
of interacting bosons evolving in the sinusoidal potential formed by a ID standing 
wave V(x) = — Vo cos 2 (fcix). In the mean field approximation, the ultracold gas can 
be represented by a single wave- function ip e (x) obeying the Gross-Pitaesvkii equation 
(GPE): 

ihMx)= f^- + V{x)+g\^{x)\ 2 ^j <k(a;) 

Projecting the solution on a basis of localised Wannier functions w„ (x), ip £ (x) = 
^2 n c n (t)w n (x) leads to the well-known tight-binding description: 

ihc n — v n c n t n —\c n —\ t n +\c n +\ -\- g \c n \ c ni (1) 

where g represents the interaction strength which is proportionnal to the s-wave 
scattering length, and we kept only nearest-neighbors hopping terms, i.e. t n+m — if 
m > 1. According to standard conventions, the hopping term is symmetric = 
t n —i= T and we measure energies in units of T and time in units of H/T. Following 
the Anderson "recipe" [T] we introduce "diagonal" disorder by randomizing the onsite 
energies v n in an interval [-W/2, W/2]; each choice of an ensemble {v n } producing a 
"realization" of the disorder, so, finally: 

ie n = v n c n Cji—i Cn+i g Cn c n . (2) 

As we use mean field theory throughout this work, we shall use the terms "interaction" 
and "nonlinearity" interchangeably. 

Spontaneous emission (SE) is inevitably - although controllable - present in an 
optical potential. For a laser-atom detuning A = ljl — 3> -To (where ujl is the 
laser frequency and ujq the atomic transition frequency and Jo the natural width of 
the transition) the amplitude Vo of the potential is Vq = hf2 2 /8A, where Q is the 
resonance Rabi frequency. The SE ratio is then given by 

_ r n 2 r n 2 

sc 4 A 2 + f2 2 /2 + r 2 /4 ~ 4A 2 

where the approximate expression corresponds to the large-detuning limit \A\ 3> 
r$,Q. For common parameters used in experimental setups, r sc < 10s _1 , T/H ~ 
1000s -1 so that in our rescaled units, the spontaneous emission rate 7 in reduced 
units is typically 7 < 10 -2 . 

We model the effect of spontaneous emission by a Monte Carlo procedure sim- 
ulating spontaneous emission events according to the probability 7. If a SE event 
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happens, as we are considering a one-dimensional lattice, we simply translate the 
momentum p of the condensate (i.e. the order parameter Fourier transform ip(p)) 
by a vector tiki, cos 9, or, in reduced units 7rcos#, with 9 being a randomly picked 
angle (in the interval [0, 2w]) representing the direction of the emitted photon with 
respect to the direction of the laser beam. Between two spontaneous emission events, 
the evolution of c„ is calculated from Eq. @. This way of modeling SE has been 
often used for individual atoms in optical lattices [18I13I1S?] and compares very well 
with experimental signals. It can however be asked if it can also be applied to the 
condensate wavefunction in the mean-field approximation. This can be justified, at 
least for low SE rates, by the following argument. In a condensate, a photon can 
be absorbed in two ways: either an individual atom aborbs it and goes to an excited 
state, or it is absorbed collectivelly generating a "delocalised" excitation. In the former 
case, the individual atom is not anymore on the condensate and does not contribute 
to the mean-field dynamics described by the Gross-Pitaevskii equation; in the latter 
case there is no changing in the condensate population and it is reasonable to con- 
sider that it is the collective mean-field wavefunction that "recoils" to compensate 
the photon momentum 0. Our model of the spontaneous emission corresponds to a 
collective recoil of the condensate wavefunction (and to a progressive increasing of 
the non-condensate fraction which is not taken into account) an approximation which 
is valid as long as this fraction stays small, that is, for low SE rates compared to the 
duration of the experiment. 

As in ref. |15) . we work with a system of fixed size of L sites, with absorbing 
boundaries to prevent wavepacket reflections. This is done by setting a slowly varying 
imaginary potential at the edges of the box. The wavepacket is thus depleted if it is (or 
becomes) large enough to touch the borders of the box, and the survival probability 
p is the integral of the wavepacket still present inside the box at time t: 



A diffusive dynamics is thus characterised by a continuous decrease of p(g, t) with t. 
In the present work typically L = 101 sites. 

In presence of interactions, the system is nonlinear and thus sensitive to the initial 
state. In order to characterise the resulting dependence of the dynamics on the size 
of the initial state, we study a particular family of initial states of square shape and 
width Lq: 



with random phases 9 n (for a discussion of the consequences of this choice, see |15|). 

3 Decoherence effects in disordered, interacting systems 

In a tight-binding approach, one associates eigenenergies and eigenstates to lattice 
sites. In a periodic lattice, each eigenenergy matches a degenerate correspondant in a 
neighbor site, which allows tunneling and diffusion. In presence of diagonal disorder, 
this degeneracy is removed, and only virtual tunneling is allowed, which leads to 
the exponentially localised Anderson eigenstates. However, spontaneous emission can 
couple Anderson eigenstates even if they are not degenrate, and thus induces diffusive 
motion in a disordered lattice. 

2 This process is somewhat analogous to what happens in the Mossbaeur effect. 



L/2 



p{g,t)= l c «W 



n=-L/2 




(3) 
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Figure 1. (Color online) Survival probability p as a function of the spontaneous emission 
rate 7 at t = 10° in the absence of interactions (g = 0). Values of the disorder W are 
2 (green triangles), 3 (red diamonds), 4 (cyan circles) , 6 (magenta stars) and 8 (yellow 
inverted triangles). The width of the initial state is Lo = 21. 

Figure [1] shows the survival probability p at time t = 10 5 as a function of the spon- 
taneous emission rate 7 for five values of the disorder W, in absence of interactions 
(g = 0). Data were averaged over a total of 500 realizations of the disorder {«„}, of 
the inital phase distribution {0 n } and of SE emission times. At time t, there have been 
on the average jt spontaneous emission events. The number of spontaneous emission 
events necessary for destroying AL depends on the disorder W. For W = 2 (green 
triangles), 10 spontaneous emissions are enough to induce complete diffusion of the 
wavepacket whereas for W = 8 (yellow inverted triangles) the wavepacket escapes 
from the box after 10 4 events (7 = 10 -1 ). 

One can interpret these results in two equivalent ways. A "static" point of view 
consists on noting that the initial state [Eq. fl3J], due to the choice of random phases, 
projects over all Anderson eigenvalues. These eigenstates have different localisation 
lengths, with a maximum at the center of the band given by £q ~ 96W~ 2 [2D], and 
zero at the borders of the band. If the disorder is small, a typical Anderson eigenstate 
has width comparable or larger than the box width L = 101, it is thus absorbed 
in a short time. The initial survival probability then corresponds to the fraction of 
Anderson eigenstates with localisation length much smaller that L. However, when a 
spontaneous emission event happens, it redistributes the quantum phases, and pro- 
jects again a part of the wavepacket on eigenstates that are either very large or 
centered in positions close to the border of the box, which are then absorbed, and 
this repeated process progressively leads to the destruction of the AL. A dynamic 
point of view considers that the initial state evolves during a characteristic localisa- 
tion time iioc ~ W~ 2 until it localises, with a typical width (a; 2 ) ~ W~ A . The phase 
redistribution due to a spontaneous emission event breaks the quantum interference 
responsible for the localisation, diffusion starts again for a time t\ oc , and then stops 
until the next SE event, and so on. For low SE ratios 7ii oc < 1, if Z? is the disordered- 
averaged classical diffusion coefficient, then one can estimate the diffusion coefficient 
induced by SE as Z?se = Dt\ oc j. 



Figure [2J presents spatial probability distributions \c n \ 2 for t = 10 5 , W — 4, Lo = 3 
for various values of the spontaneous emission rate 7. For 7 = (blue solid line) this 
distribution is exponentially localised |c„| 2 oc exp(— 2 |n| /£), where I is the localisa- 
tion length, comparable to £q(W), the maximum localisation length of the Anderson 
eigenstates. For 7 = 10~ 5 (green dotted line), the system is only weakly affected by 
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Figure 2. (Color online) Density distribution \c n \ (t = 10 ) vs site index n, for g = 0, W = 
4, Lq = 3. Values of 7 are : (blue solid line), 10~ 5 ( green dotted line), 10 4 (red dashed line) 
and 10 -3 (dot-dashed cyan line). Anderson localisation, characterised by the exponential 
(triangular in logarithmic scale) shape is progressively destroyed when decoherence increases. 

spontaneous emission (as just one event happens in average during the time evol- 
ution), a small diffusion is observed in the tails of the distribution. For 7 = 10 -4 
(red dashed line) and 7 = 10~ 3 (dot-dashed cyan line), the distribution has taken a 
Gaussian shape, indicating that the decoherence-induced diffusion has become dom- 
inant. Alternatively, this process can be interpreted "statically" as a phase reshuffling 
provoked by spontaneous emission, which excites eigenstates centered everywhere in 
the box L, so that the localisation length is not anymore a relevant length scale in 
the problem. 



Consider now the repulsive interacting case (g > 0) . It is useful, if not completely 
rigorous, to interpret the nonlinear term in Eq. ([%]). as a dynamical perturbation 
w NL = g\c n \ 2 of the on-site energy v n , depending on the population of the site. For 
v nl ^ _ which i m phes a low-density, spatially extended, wavepacket - Anderson 
localisation is expected to survive for a very long time: Anderson eigenstates whose 
localisation length i is inferior to the size of the box remain in the system and give 
non-zero survival probability. For w NL ~ W, if the c„(i) vary rapidly, the dynamical 
correction can bring temporarily neighbour levels close to degeneracy, restablishing 
transport. This mechanism is most efficient in the so-called chaotic regime, where the 
c n (t) display a stochastic character. If v Nh ^> W, neighbour sites are decoupled from 
each other even for small population differences, which leads to a different, nonlinear, 
type of localisation called self-trapping, favored for narrow (small Lq) initial states. 
Figure [3] displays the survival probability p(t = 10 5 ) as a function of the interacting 
strength g, for W = 4 and four values of the decoherence level 7 from to 10 -3 . 
The three dynamical regimes are clearly visible, a plateau for small values of g due 
to surviving AL, a dip created by the chaotic diffusion for intermediate values of g, 
and a reduction of diffusion due to self-trapping for high values of g, leading to a 
complete freezing for g w 300. Clearly, the three dynamical regimes are not affected 
in the same way by decoherence: Anderson localisation, which relies on delicate phase 
relations is more sensitive to decoherence. Comparison of the curves for 7 = (solid 
blue) and 7 = 10~ 5 (dotted green) shows that even one SE event in the average has 
an observable effect on AL. In the chaotic regime, the diffusion due to chaotic motion 
and diffusion due to decoherence tend to add to each other, but a strong diffusion due 
to decoherence decreases the site populations, and thus the nonlinear term gJ2 n l c «| 2 - 
This effect of dilution due to decoherence explains the fact that the frontier between 
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Figure 3. (Color online) Survival probability p as a function of the interacting strength g 
for t = 10 s . Values of 7: (blue solid line), 10~ 5 ( green dotted line), 10 4 (red dashed line) 
and 7 = 10 -3 (dot-dashed cyan line). Three dynamical regimes are clearly visible: Localized 
(low g), chaotic (intermediate values of g), and self-trapped (high g). Other parameters are 
W = 4 and L = 3. 



the various regimes is shifted to higher values of g when 7 increases. Self-trapping, not 
relying on quantum inteference, but only on populations |c n | 2 (which are not directly 
affected by SE) tends to be insensitive to decoherence: for g « 320 all curves converge 
to unitary survival probability. 

Figure 0] shows spatial distributions at t = 10 5 for g = 10 (a) and g = 320 
(b). In the chaotic regime (a) and for low decoherence rate (solid blue, dotted green 
lines), the wavepacket has its exponential shape preserved in the center of the box 
but non-exponential tails appear. This is due to the fact that interactions depopulate 
the center of the box (where populations and thus nonlinearity is maximum), but it 
preserves the exponential shape |15j . The question whether Anderson localisation is 
destroyed in this case is tricky: The wavepacket spreads along the box but preserves 
its typical exponential-shape. Adding decoherence to the system gives a far more solid 
answer to this question. One sees in Fig. 2] that the exponential tail peak is almost 
completely destroyed after 10 spontaneous emissions on the average (dashed red and 
cyan dot-dashed lines) , which strongly suggested that it is indeed due to interference 
effects and can thus be attributed to a survival of AL. On the other hand, as expected, 
the self-trapping regime Fig. 0b is very robust against decoherence. 



4 Scaling laws 

Sensitivity to the initial conditions is a common feature of nonlinear systems. The 
nonlinear correction v NL obviously depends on the local density and is therefore highly 
sensitive to the initial width of the packet: the linear regime is favored for high values 
of Lq whereas the chaotic regime and self-trapping regime are expected for spatially 
concentrated wavepackets. In ref. |15) . we showed the existence of scaling laws on the 
width Lq of the initial wavepacket, which in particular allowed us to classify the three 
regimes as a function of a scaled interaction strength g = gL^ 11 with s = 0.76 ± 0.08. 
Figure [5^, compares the survival function p(t = 10 5 ) for different values of the initial 
state width Lq; one clearly sees that the crossovers between the quasilocalised, chaotic 
and self-trapped regimes indeed strongly depend on L - Figure [5b displays the same 
curves plotted in terms of g, showing that the crossovers become independent of the 



Will be inserted by the editor 



7 




Figure 4. (Color online) Spatial distributions at t = 10 5 for g = 10(a) and g = 320(b), 
and 7 = (solid blue line), 7 = 10" 5 ( green dotted), 7 = 10 4 (red dashed) and 7 = 10 3 
(dot-dashed cyan). For low nonlinearities and low decoherence levels, Anderson localisation 
survives for long times, but is finally destroyed by decoherence. The self-trapping regime is 
almost insensitive to decoherence. Other parameters are W — 4, Lo = 3. 
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Figure 5. (Color online) Survival probability as a function of the nonlinearity for W = 1, 
7 = 10 _J . (a) No scaling, (b) The scaling in g — gL^ s (s = 0.76) allows a clear identification 
of regions corresponding to localised (g < 0.1), chaotic (0.1 < g < 10) and self-trapped 
(g > 10) behaviour. Values of Lo : 3 (blue squares ), 7 (green triangles), 13 (red diamonds), 
21 (cyan stars), 31 (magenta circles), 41 (yellow inverted triangles). 

initial state, even in the presence of spontaneous emission. As opposed to interactions, 
decoherence effects are not sensitive to the initial state so it is not surprising that the 
renormalized nonlinearity remains a "good" parameter. 

Figure [5] shows that the probability p(t) can also be approximately scaled in Lq, 
but the scaling depends on the dynamical regime: p = pL$ with v = 0.52 in the 
chaotic regime and v = 0.31 in the self-trapped regime (and, trivially, v = in the 
quasilocalised regime) . 

5 Conclusion 

We have thus characterised the influence of decoherence on the dynamics of interacting 
bosons in a disordered potential. Spontaneous emission changes considerably the time 
behaviour of the survival probability, but the use of the scaled parameter g takes into 
account the crucial sensitivity of the system to the initial conditions, even in presence 
of decoherence induced by spontaneous emission, allowing a global characterisation 
of the dynamics. This shows that the effects of interactions and decoherence tend to 
add up rather than compete. These above results put into evidence the usefullness of 
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Figure 6. (Color online) Scaling of the survival probability p(t — 10 5 ) (log- log plot) as a 
function of the scaled nonlinearity g. (a-c) 7 = 10 -5 and W — 4, (d-f) 7 = 10 -3 , W = 1. 
Left column: p not scaled and full range in g; middle column: p scaled in L®' 52 zoomed on 
the chaotic region 0.1 < g < 10; right column: p scaled in L[j' 31 zoomed on the self-trapping 
region g > 10. Values of Lo : 7 (green triangles), 13 (red diamonds), 21 (cyan stars), 31 
(magenta circles), 41 (yellow inverted triangles). 



the scaling 011 the initial state as a tool allowing state independent characterisation 
of the dynamics, which thus constitutes an important element of a "new language" 
better adapted to describe the quantum mechanics of nonlinear systems. 

Laboratoire de Physique des Lasers, Atomes et Molecules is Unite Mixte de 
Recherche 8523 of CNRS. Work partially supported by Agence Nationale de la Recher- 
che's LAKRIDI grant and "Labex" GEMPI. 
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